Analytical kinetics of clustering processes with cooperative 
action of aggregation and fragmentation 
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Some models of clustering processes are formulated and analytically solved employing generating 
functions methods. Those models include events which result from combined action of the coag- 
ulation and fragmentation processes. Fragmentation processes of two kinds, so-called similar and 
arbitrary, ones, are brought forward, and the explicit forms of their solutions are produced. This im- 
plies some possibility of existence of different aggregation mechanisms for clusters creation differing 
in their inner structure. All the models are based on "the three- level bunch" scheme of interaction 
between the system states. Those states are described in terms of the probability to find the system 
in the state with an exactly given number of clusters. The models are linear in the probability 
functions due to the assumption that the rates of elementary acts are permanent. Some peculiari- 
ties of application of the generating function method to solution of the linear differential-difference 
equations are revealed. The illustration of the problem in terms of a traffic jam picture is not a 
specific one. 
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I. INTRODUCTION 



In the present work we develop and define more accu- 
rately our earlier considerations A cluster is gener- 
ally understood as either a number of things of the same 
kind growing together or number of particles, objects, 
etc., in a small, close group. This idea is a very gen- 
eral and not a specified one. An aggregation process in 
the form au + bu = (a + b)u is referred as coagulation, 
when coefficients a, b represent the quantities of the scale 
unit of u which are coalesced with time. Aggregation 
and fragmentation are a couple of mutually inverse pro- 
cesses. The physical scales ( spatial, temporal, value of 
masses, e.t.c.) vary in the many orders of a magnitude 
for such a processes. That is why, an idea of existence 
of a universal description of the above phenomena arises 
and the unification of both direct and inverse processes 
in a general class of clustering processes seems being of 
a natural one. 

The clustering processes resemble to a certain degree 
the nucleation processes. This roots in the mathematical 
description being general for kinetics of such processes 
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despite the lack of obvious resemblance in their actual 
more precise details. On the one hand, we can see those 
approaches in the course of investigations on the molecu- 
lar and submolecular level, in theories of condensed mat- 
ter, nuclei and nuclear chains |2j. On the other hand, 
clustering of disperse systems are considered in astro- 
physics (forming of cosmic objects), atmospheric science, 
chemistry, ... j3j. 

For example, an expanding universe is formed not at 
once. Clusters grow by coalescence of smaller clusters. 
Their growth kinetics is like to the kinetics of coagula- 
tion. In what follows we formulate basic equations and 
outline the methods for their solution. Moreover, one 
could expect that theoretical tools developed to describe 
physical systems can be exploited in other fields, such as 
ecology of computation Q] or biology, economics, trans- 
port problem, etc. & |(|. 

Consider the kinetics of formation of a G-cluster using 
the picture of a one-way motor lane. We assume that 
the starting configuration is G independent cars on the 
motor lane, the leading one being the slowest, and no one 
can pass over each other. So, each initial cluster contains 
one car. 

The process begins at t = 0. On passing some time i, 
the initial G cars group in clusters containing gi , g% . . . g s 
cars. These clusters go on to coalesce. The problem 
is to determine the time evolution of the probability 
w s (gi 1 g2 ...g s ',t) to find s clusters g\, gi . . . g 8 . The sum 
of their masses gk are subjected to the constraint (con- 
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servation law): 



E 9k 
fe=i 



G 



(1) 



in the system. Thus, we concern only of non-relativistic 
events and study systems with a permanent (additive) 
mass and a finite number of particles. 

Our goal is to formulate and investigate the exactly 
solved models of the clustering (dissociation) processes, 
including those which result from combined action of a 
certain aggregation and fragmentation. 

The paper is organized as follows. In Section 0, we 
explain the stochastic motion of our objects, obtain the 
probability of finding out the system in the state of ex- 
actly s clusters and dependent on the time average num- 
ber of clusters by means of introducing a generating func- 
tion. Later on, we formulate a master equation governing 
the time evolution of the probability of finding the clus- 
ters of various masses (Sect. |[Il|) . We solve that problem 
applying the Laplace transformation with respect to t to 
master equation (Sect. IV). Then we find the probabil- 
ity to detect a cluster of assigned mass g by summation 
of w s over all g k irrespective of the distribution of other 
participants, except for the selected one (Sect.0). Some 
properties of those convolutions are true due to the iso- 
morphism between a set of generating functions with a 
product operation and a set of w(g) with convolution (see 
Appendix as well), which have the structures of semi- 
groups with unit. In the following we consider processes 
of similar fragmentation (Sect. VI), combined act ion of 
aggregation and similar fragmentation (Sect. VII ), an d 
the process with an arbitrary fragmentation (Sect. VIII). 
In conclusion, we discuss our results. 



In particular, if initially there were exactly G indepen- 
dent cars, the function Wq(s) obeys the equation 



Wq(s) = A(s - G), 



(4) 



with A being the Kronecker delta: A(0) = 1, and A = 
otherwise. 

Equation (^) can be solved by introducing the gener- 
ating function: 



F(z,t) = J2W(s,t) 



y S— 1 



Combining Eqs. (||) and @ gives: 

d t F = (1 - z)d z F. 



(5) 



(6) 



The rate 7 is included into the definition of time. The 
initial condition for an initially monodisperse system is 
rewritten in terms of z as: 



F(z,0) = z 



G-l 



(7) 



The solution of Eq. (|6|) with the initial condition, Eq. 
(0), has the form: 



F(z,<) = [l-e- 4 (l-z 



1(3-1 



(8) 



The probability W(s, t) is thus expressed in terms of bi- 
nomial distributions: 



W(s, t) = C^e-^-^il - e-*) 



-t\G-s 



(9) 



It is no problem to find the time dependence of the av- 
erage number of clusters: 



a{t)=d z zF(z,t)\ z=1 = l + (G-l)e- 



(10) 



II. NUMBER OF CLUSTERS: PURE 
AGGREGATION PROCESS 

Let 7 be the rate of an elementary coalescence act; 
two adjacent clusters produce a single one (for instance, 
a dimer is formed when a car catches up another one). 
We assume that 7 is ^-independent. Then we can char- 
acterize the situation by the number of intervals between 
adjacent clusters. If there are s clusters in the system, 
the number of intervals is s — 1. Each coalescence act 
annihilates one interval. The number of ways to do this 
is exactly equal to the number of intervals. 

Let W(s, t) be the probability of meeting exactly s 
clusters at the time t. Then: 

M^=rts W ( il+1) t)-(s-l)W(8 ) t)]. (2) 

One can observe a "three- level bunch" scheme of transi- 
tions between the three nearest (on s) states of the sys- 
tem under consideration. Equation (Q) should be supple- 
mented with the initial conditions: 

W(s,0)=W (s). (3) 



III. MASS DISTRIBUTION IN A PURE 
AGGREGATION PROCESS 

In analogy with the kinetics of disperse systems, we 
shall refer to gu as the cluster mass. Our goal now is to 
formulate the master equation governing the time evolu- 
tion of the probability w s (<?i , g% . . . g s ; t) to find the clus- 
ters of masses 31,32 ■■ ■ at the time t. This equation is 
formulated as follows: 

J*jg»= £ w a+1 (gi...g' k ,g' k+1) g' k+2 ...g' s+1 ) 

[9] ' k xA(g k -g' k -g' k+1 ) ( n ) 
x A (5fc+2 ~ 9k+i) ■ ■ ■ Mg's+i - 9s) - (s - l)w s . 

The meaning of the terms in the r.h.s. of Eq. (|ll| ) is 
rather apparent. The rate of losses is simply propor- 
tional to the number of empty intervals (the rate con- 
stant 7 is included into the definition of time). The gain 
occurs each time when two clusters of masses g' k and g' k+1 
coalesce producing a new cluster of mass gu- Other A-s 
simply restore the serial numbers of gi clusters with i < k 
for the system of s clusters. 
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Of course, initial conditions to Eq. should also be 
specified. We again assume that initially there were G 
separate cars: 

«te(l,l...M = 0) = l, (12) 
and all other probabilities are 0. 

IV. SOLUTION TO BASIC EQUATION 

On applying the Laplace transformation with respect 
to t gives, instead of Eq. (Thy), 



(p + s 



l)w a (gi,92 ■ ■ ■ ) 
• 9k>9k+n9k+2 ■ 



■9s + i) 



E 

xA(g k -g' k -g' k+l ) 

xA (g'k+2 - 9k+i) ■ ■ ■ Ha's+i - 9s), 



(13) 



where barred w stands for the Laplace transform of 
w(gi, 32 ■ • ■ ',t). The last equation of this set is readily 
solved (Eq. (ft)) to yield: 



w G 



P+G-l- 



(14) 



Now let us try to seek for a solution to Eq. (|llj) in the 
form: 



w a (gi,g 2 -- - ;,p) 



A B (g 1 ,g 2 ... ) 

(p+G-l)(p+G-2)...0+s-l) ' 



(15) 



where the coefficients A are independent of p and satisfy 
the following set of recurrence relations: 



A s (gi,g 2 ...) 
J2 A s+1 { gi . . . g' k , g' k+1 , g' k+2 
la']* 

xA( 5fc - g' k - g' k+1 ) 
xA (9k+2 ~ 9k+i) ■ ■ ■ Mg's+i ~ 9s) 



g's+i) 



(16) 



A useful sum rule 



originate from the interchange of generations at a given 
G, and the proper weights depend on the whole set of 
possible transitions from s + 1-states to the s-state under 
consideration. 

With the induction method, one obtains that: 

• Ag(1, ■■■ , 1) = 1 is the only possible value; 

• the number of terms in formula Eq. ( |l6|) is equal 
to 



E 

k=l 



9k - 1 = 9k - ^ 1 = G - s 

k=l k=l 



for each s fixed and k fixed. Under the 
inductive assumption, one can write down 
A s+ i(g[, . . . ,g' s+ i) = A s +i. From this it follows 
that A s (gi , . . . ,g s ) = (G — s)A s+ i irrespective of a 
specific distribution of g±, . . . , g s . 

The recurrence equations obtained just now 



have solutions 



As = (G-s)A s+1 , 
Aq = 1 



A s = (G- s) 



(19) 



(20) 



The time dependence can be readily restored by using 
the inversion 

l 



(p+s-l)(p+s)...(p+G-l) 

♦^e-^Cl-e-*) -- 



(21) 



Equations (16), (|l7|) , and (|2ll ) reproduce Eq. @ as well. 
The final result is formulated as follows: 



w s (gi,g2, 



;i) = e-( s - 1 ) t (l-e~*) 



-t\G-s 



xA(G-gi - 32 



(22) 



V. SINGLE-CLUSTER DISTRIBUTION IN A 
PURE AGGREGATION PROCESS 



Q s = sQs+i, 



(17) 



follows immediately from Eq. (|16|), where Q s = 
^2A s (g\,g2--- g s ) (summation runs over all g), or 



(G-l)' 



(18) 



In fact, the expression A s (gi,... ,g s ) depends on s 
only. It does not depend on the distribution of num- 
bers gi, . . . ,g s provided that J2k=i 9 k — G, gk > 1 are 
conserved. 

It can be seen from Eqs. @, (||), (|), and (|) that 
the problem under consideration splits into two subprob- 
lems. The first one is the time evolution problem. It 
deals with transitions between different states of the sys- 
tem and connects each three nearest adjacent states. The 
second subproblem is to scrutinize mass spectra. It is a 
pure combinatorial task. In fact, we have to deal with 
some population dynamics. Mass spectra at instant t 



Here we determine the probability to find a cluster of 
mass g irrespective of the distribution of other partici- 
pants. To this end we sum w s over all g k except one (<?i, 
for example): 

w{g,t) =T, w s(g,92, ■ -.gait) 



= e -( S -l)t( 1 _ e -i)G- s 
x E MG-. g- : g 2 --..gs)- 



Using the identities 



A(?) = 



1, 9 = 
0, q=l,2,. 



JL f J*L 
2-iri J z*z' 



0, q = 

1, 9= 1,2, 



1 r zdz 



2m J z«(l— z) ' 
1 X dz (~ir /-iR 



(23) 



(24) 
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A (6-) 



FIG. 1: Generation scheme for G=6 



one finds the convolution in Eq. (£2 



-t\G-s 



w(g,t) = e- ( - s -V t (l-e- t ) 

v J— £ z"- 1 dz 

X 27TI T ZV-S(1- Z )°- 1 Z 

= C S ~ 2 ^-( S - 1 )*C1 _ -t\G-s 



(25) 



, G _ g _ ie -(°-v\i-a-r 

Some important properties of the convolutions are dis- 
cussed in Appendix. 




FIG. 2: W-functions versus time for pure aggregation process. 
G = 5, 7 = 1. All curves W(s, t) are marked by the number a 
of the proper state. 



VI. SIMILAR FRAGMENTATION PROCESS 

Let us consider a process of pure fragmentation (disso- 
ciation, decay) of clusters. We assume the inner G-cluster 
structure at t — to be similar to the picture at t = 
with G cars of unit mass in Sect. [0] and, thus, with G — 1 
intervals between them . An analogous assumption ap- 
plies to all the s-clusters at t ^ 0. We will understand 



" the similarity of an inner cluster structure and an outer 
one" in such a purport. In addition, let us suppose a 
qualitative equivalence of all the inner intervals between 
constituents of a cluster. 

At the above pure aggregation mechanism the place 
where a coalescence act has happened is partly forgotten. 
When two adjacent clusters enumerated k s and (k + l) s 
according to enumeration of an s-cluster state with sizes 
(see below) (g'k)s and (g*k+i)s coagulate, a new cluster 
of the size (gk)s~i = (g , k)s + (<?'fc+i)s arises. The cluster 
ordinal number fc s _i is given in terms of the new (s — 1)- 
cluster state which originates from the act of coagulation 
as in Sect. III. 



Evidently, a partial loss of memory on the way in which 
the microscopic state has been created comes about at 
the very place because the inverse problem of one-to- 
one rebuilding of the previous s-cluster state (the outer 
state) cannot be solved. On the other hand, the things 
are the same with the interior structure of some cluster. 
Its cluster parents cannot be reconstructed one-to-one as 
well. When one looks at a system state or some cluster 
as something given, the exact information about ordinal 
numbers and sizes of adjacent cluster parents and even 
the mother state is forgotten. That is why one can talk 
about a loss of memory in such a process. In the whole 
it is the reason to postulate some equivalence of con- 
stituents and intervals between them in the inner cluster 
space. 

This set of assumptions results in the dependence of 
the probability functions on the size of a cluster and the 
time only. 

Let us realize the cluster size as the number of particles 
confined into the given cluster, let 7 be the rate of an 
elementary fragmentation act. 

If a fragmentation rate is proportional to the clus- 
ter size minus unit, i.e., the number of possible rupture 
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places is equal to the number of inner intervals, the equa- 
tion 

= 7 [(G - a + - 1, t) - (G - s)W(s, t)}, 

(26) 

with initial conditions 

W(1,0) = 1; W(s,Q) = 0, ifa^l, (27) 

describes the process under consideration. The r.h.s. of 
Eq. ( |2li| ) consists of a gain term due to decay of clusters 
belonging to an s — 1-cluster state and a loss term due 
to decays of those clusters belonging to s-cluster state 
that produces the clusters pertaining to an s + 1-cluster 
state. Equation ( p6| ) can be solved by using the generat- 
ing function introduced by the equation 



&F(z, ()-h 2 7 + ^7)|F(z, f) 
= j(zG- z-G+l)F(z, t), 

with the solution 

-z-el-i'l+el-T') 2 ' 



F(z>*) = -- 



(28) 



(29) 



Of course, one recognizes a usual Poissonian process here. 
It seems contextual and, hence, quite reasonable to name 
such a process as a similar fragmentation (process). 
For example, for G = 5 



Wi 



3 (-4 7 i) 



W 2 =4e<- 3 T*) -4e<- 4 7«), 



(30) 



(31) 



W 3 = 6e(" 2 ^) -12e(- 3 ^)+6e(- 4 ^), (32) 



with initial conditions 



W(G, 0) = 1; W(s, 0) = 0, if s + G. 



(37) 



The r.h.s. of Eq. ( pq ) consists of gain terms due to coag- 
ulation of clusters from an (s + l)-cluster state and dis- 
sociation of those clusters belonging to an (s — l)-cluster 
state and loss terms due to simultaneous coalescence and 
dissociation of clusters belonging to an s-cluster state. 
To make things more clear, we could rewrite Eq. ( j36|) in 
the form 

dW(s,t) 

= { llS W{s + l,t)+^(G - s + l)W(s - l,t)) (38) 
-(7i (« ~ l)W(s, t) + 72 (G - s)W(s, t)). 

These equations can be solved by using the generating 
function F(z,t) defined by the equation 



f t F(z,t) 

-(-z 2 7 2 -271 +z 72 -l-7i)^F(z, t) 
= l2 {zG-z-G+l)F(z, t), 



whose solution is 



F(z,t) 



' Z 72+71 e< ' 7| z-7ie' 1T '+7l ^G / 

(272 + 71 e ( ~' 7) z -71 e ( ~* 7) + 71), 



where 

7 = 71 + 72- 
For example, for G = 5 



(39) 



(40) 



(41) 



W 4 - 4e^ ( ) - 12e(- 2 ^) + 12e^^ - 4e(" 4 ^), 



(33) 



Wi = .6830134554 (e^ 11 *) - l) 4 , 



(42) 



W 5 = l-4e(-T*) +6e(- 2lt ) -4e(" 3 7*) + e (- 4 ^). 

(34) 

The average number of clusters reads 

g(t) = d z zF(z,t)\ z=1 = e (-T") + G(1 -e^*)). (35) 



W 2 = -2.732053822 (e^ 1 - 1 *) - l) 3 (.1 + e^ 1 - 1 *)), 



(43) 



W 3 = 4.098080732 (e^ 1 - 1 *) - l) 2 (.1 + e^ 1 - 1 *)) 2 , 



(44) 



VII. PROCESS OF AGGREGATION AND 
SIMILAR FRAGMENTATION 

Let us consider such a clustering process, which runs as 
a result of some combined action both of the aggregation 
and the similar fragmentation. Let "/i and 72 be constant 
rates of an elementary coalescence act and an elementary 
fragmentation act, respectively, 

dW(s,t) 

dt , 

= n[sW(s + l,t)-(s-l)W(s,t)] (36) 
- 72 [(G - s)W(s, t)-{G-s + l)W(s - 1, t)], 



W 4 = -2.732053822 (e^ 1 - 1 *) - 1) (.1 + e^" 1 - 1 *)) 3 , 

(45) 

W 5 = .6830134554 (.l + e^ 11 *)) 4 . (46) 
The average number of clusters is given by the formula 
5(f) = d z zF{z,t)\ z=1 

= 1 _ 72+71 , G (72+71 e'-' 1 ') i 47 ) 
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It should be noted, that the mass spectra of clusters for 



FIG. 3: W-functions. Process of aggregation and similar frag- 
mentation. G = 5, 7 = 1.5,7i = 1,72 = 0.5. All curves 
W(s,t) are marked by the number s of the proper state. 



process Eqs. ( [48[ ) and ( |49| ) differ from those in Sect. Ill 
r-K 

Equation fl48|) can be solved for each G. If, for example 
G=5, 



W{s,t) = Y:% 1 C i e^ t A sU 

C = [.44721, -.44721, -.44721, -.44721, -.44721], 
A = [-10.95389, -5.73117, -2.57163, -.74329, 0], 



.01096, .11376, -.43096, .77537, -.44721 
-.10909, -.53825, .67732, .19904, -.44721 

.42839 , .67814 , .36055 , -.16309 , -.44721 
-.77746 , .19355 , -.15969 , -.36411 , -.44721 
.44721 , -.44721 , -.44721 , -.44721 , -.44721 



(50) 

(51) 
(52) 

(53) 



VIII. PROCESS WITH AN ARBITRARY 
FRAGMENTATION 

Let us assume that clusters lose the memory of their 
aggregation history completely. This means that a hypo- 
thetical mechanism of such a kind aggregation loses not 
only the place of a coalescence act (Sect. VI), but in this 
case the inner structure of a cluster does not coincide 
with outer one, becomes a uniform one, and let us say 
"is beside itself. Shortly speaking, the mass distribu- 
tion inside a cluster becomes a uniform one under such 
an aggregation, i.e., the inner structure does not depend 
on the cluster size (see Sect. VI, a nd h ence, the actual 



number of inner intervals (see Sect. VII) is equal to 0. 

However, let us apply the "three- level bunch" scheme 
of transitions (connections) between the three nearest 
cluster-states, widely used above. That is, in fact, one 
of the main subjects under consideration in this paper. 

In this case, it is reasonable to assume that the frag- 
mentation rate depends solely on the outer cluster state 
indices. I.e., the rate depends on the integer number of 
the clusters or the intervals minus unit between them for 
the states belonging to such a three-level bunch. Adding 
the contextual "natural" condition of stopping the frag- 
mentation process at s = G, the corresponding equation 
reads 



= 7l [sW(s + l,t)-(s- l)W(s, t)} - 72 
x [(1 - A(G - s))sW{s, t)-(s- l)W(s -l,t)}, 



with initial conditions 



W(G,0) = 1; W(s,0) = 0, ifs^G, 



(48) 



(49) 



where 7 i and 72 are the constant rates of an elementary 
coalescence act and of an elementary fragmentation act, 
respectively. Equation ( |48l ) is constructed so as to absorb 
the conditions W(G + 1, t > 0) = 0, W^t > 0) = 0. We 
shall see below, the solutions of Eq. (j48|) bear out this 
assertion. 




FIG. 4: W-functions. Aggregation with an arbitrary frag- 
mentation. G = 5, 7 = 2, 71 = 1, 72 = 1. All curves W(s,t) 
are marked by the number s of the proper state. 



IX. CONCLUSION 

In this paper, we are interested in kinetics of cluster 
formation and construct a number of analytically solved 
models of such a process. It does not matter whether a 
motor lane, a computational network, space dust or glass 
is a real environment to match such a script. The illus- 
tration of the problem in terms of a traffic jam is not a 
specific one. The known technical and natural phenom- 
ena of traffic jam or aggregation (coagulation) and their 
evolution have been considered as processes of clustering 
(nucleation) . 

Brought forward in this work, the kinetics models take 
account of both the aggregation (coagulation) and the 
fragmentation (decay) processes. Previous studies of 
nonlinear Lotka-Volterra systems Q| brought us to a 
search for a possibility to give a linear description of 
those very complicated and nonlinear situations or for 
something akin to such a picture. A dynamical descrip- 
tion of some system could be substituted by a stochastic 
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one. We consider one-dimensional cases (in the coordi- 
nate sense) but the scenarios should be valid for the three 
dimensions if there is a spherical symmetry. 

When we deal with a system of a finite number of parti- 
cles, a natural way of the above substitution is to use the 
language of enumerated states of that system. A state of 
that sort is characterized by the population number and 
probability function to reveal the system itself in this 
state exactly. The "three- level bunch" scheme of tran- 
sitions between the three nearest states of the system is 
natural as well. Of course, that probability should de- 
pend upon the probability of something else to happen. 
In the case considered, it could be an act of coalescence 
or decay ( dissociation, fragmentation, etc., where the 
term depends on the application ). Thus, we have just 
met a product of probabilities. What can one do? 

The probability (rate) of the above acts could be de- 
pendent or independent on the system states or its partic- 
ular attributes. In the case of such a dependence, one can 
say nothing without a special investigation. On the con- 
trary, the rate independence of the above circumstances 
makes the situation a linear one in the state probability 
function. Let us take into account that the rates may 
be independent not only of macro parameters but also 
of the cluster size (e.g., if coalescence/fragmentation de- 
pends only on the valency of some chemical clusters). 
This leads to our key mathematical assumption of the 
above rates being constant. That is why we obtain linear 
analytically solved master equations of that clustering 
kinetics, which are also evolutionary- type equations. 

The very method to solve most of those problems is in 
use of the generating function. 

But some peculiarities of applying of this method for 
solving the linear differential-difference equations are re- 
vealed. It becomes useless when the structure homogene- 
ity of terms in the r.h.s. of equations discussed is violated 
as a result of combining aggregation and fragmentation 
terms, i.e., owing to violation of the Markovian semi- 
group structure of those r.h.s. The latter, as well as 
symmetry, momenta and other questions (e.g. univer- 
sality of the three-level bunch scheme ) not included in 
this sketch will be considered in more detail in papers 
to follow, as well as some details of the utilization of a 
MAPLE- program, used to obtain some analytical results, 
which is to be published in || . 

Other questions about the mass spectra and the single- 
cluster mass distribution in a pure aggregation process 
help to display the combinatorial nature of the problem 
under consideration and its closeness to some popula- 
tion dynamics problem. They are solved by means of the 
Laplace transformation, which leads to the calculation 
of certain convolutions. The properties of those convo- 
lutions are true due to the isomorphism between a set 
of generating functions with a product operation and a 
set of w(g) with a convolutions, which results in a useful 
formula. 

It should be remembered, that we considered two types 
of fragmentation: the similar fragmentation and the ar- 



bitrary one. They lead to different results. Such a dif- 
ference can be realized as a manifestation of distinction 
between collective and more discrete additive properties 
of a cluster on a very abstract level, far off a specific na- 
ture of that cluster. Moreover, we shall see the possible 
existence of several aggregation mechanisms (for a variety 
of reasons). Thus, one and the same stochastic descrip- 
tion of the aggregation can mask a number of processes 
diverse in their dynamical details. 

However, the generality of the results is restricted. The 
conjecture that the rates are constant (71,2 = consti^) 
is nowadays somewhat at odds with what is known 
about molecular and atomic clusters, and processes for 
nanoscale units of matter. Obviously, the mentioned con- 
jecture reduces a mathematical universality as well. We 
are grateful to the referee for drawing our attention to 
the necessity of making this remark. 

Nevertheless, the authors hope that appropriate pro- 
cesses of such a type may be found by the following 
reason. The presented formalism is independent of any 
physical scale. The only conjecture of the possibility to 
draw distinction between certain particles and intervals 
between them has been done. Such an admissibility could 
be combined with not only a classical conception but a 
quantum one under some conditions (e.g., if the first Born 
approximation is valid ) as well. 

Present work belongs to the stream originated from 
the famous Smoluchovski's articles S. 
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APPENDIX A 

Let 

JV = 0,1,2,..., 
a(n) = a n , n e A^, 

(Al) 

A = {a : N Q -> R}, v ' 

F = {f:C^C,f(z)= £ a n z n }. 

n£N 

The reversible mapping Z : A — > F 

f = Z(a) so as f(z) = a{n)z n (A2 \ 

neN V ; 

has the inverse mapping Z~ x 

a = Z-\f) so as a(n) = § f (z) ^ , 

(A3) 

n 6 N . 
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The convolution o is an operation on a set A such that the unit A: 

k 

c = aob iff c(k) =J2a(k- j)b(j). (A4) 

3=0 



A(0) = 1, A(n) =0, n + 0, 

(A6) 

Z(A) = 1, Z-^l) = A. 



The mapping Z is a morphism of a semigroup (A, o) 
to a semigroup (F, •). Z maps the convolution o to the Here is a useful formula 
product • : 

(di o a 2 o . . . o a n )(mj 

c = ao()^Z(c) = Z(a)Z(!)) (A5) = E a i(h) ■ ■ ■ a n (j n )A(m - ji - . . . - j n ) (A7) 



The associative and commutative semigroup (A, o) has 
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